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In a recent experiment, Lai et al. [Phys. Rev. B 75, 033314 (2007)] studied the apparent metal- 
insulator transition (MIT) of a Si quantum well structure tuning the charge carrier concentration 
n. They observed linear temperature dependences of the conductivity a(T, n) around the Fermi 
temperature and found that the corresponding T — > extrapolation ao(n) exhibits a sharp bend 
just at the MIT. Here, reconsidering the data published by Lai et al, it is shown that this sharp 
bend is related to a peculiarity of a(T — const., n) clearly detectable in the whole T range up to 4 K, 
the highest measuring temperature in that work. Since this peculiarity seems not to be smoothed 
out with increasing T it may indicate a sharp continuous phase transition between the regions of 
apparent metallic and activated conduction to be present at finite temperature. Hints from the 
literature of such a behavior are discussed. Finally, a scaling analysis illuminates similarities to 
previous experiments and provides understanding of the shape of the peculiarity and of sharp peaks 
found in dlog 10 a/dn (n). 

PACS numbers: 71.30.+h,73.20.Fz,73.40.Qv 



I. INTRODUCTION 

The question whether electronic conduction in two- 
dimensional disordered systems is exclusively nonmetal- 
lic, so that the resistivity p(T) always diverges as temper- 
ature T — > 0, or whether it can also have metallic charac- 
ter has been under controversial debate for three decades. 
The existence of a corresponding metal-insulator transi- 
tion (MIT) was denied by the localization theory of Abra- 
hams et al. which however neglects electron-electron in- 
teraction. Thus it came as a surprise when Kravchcnko 
et al., who had varied the charge carrier concentration 
n in high mobility MOSFET samples, first reported a 
strong decrease of p down to 20 mK^ They considered 
the conduction in the respective (T, n) region as metal- 
lic, an interpretation which was questioned in particular 
by Altshuler and Maslov^ compare also Rcfs. 1 andU 
The nature of the apparent metallic state has remained 
puzzling up to now, for recent reviews see Rcfs. 

Provided the apparent MIT is a real phase transition at 
T = 0, it can be considered as a quantum critical point 
Commonly, the transition is expected to be smoothed out 
at finite T. However, this smoothing is not a must for a 
quantum phase transition^ For example, the Ising fer- 
romagnct LiHoF4 in transverse magnetic field undergoes 
a quantum phase transition at T = 0, where the quan- 
tum critical point is the beginning of a line of continuous 
phase transitions at finite Thus both the situations 
sketched in Fig. 1 are conceivable. It is the aim of this 
work to uncover indications in the literature for the MIT 
in two-dimensional systems being likewise connected with 
a line of continuous phase transitions at finite T. 

The present study is motivated by the scaling of the T 
dependences of p for various n, p(T,n) = p(T/To(n)), 
reported by Kravchenko et al. for MOSFET's.^ For 
n < n c , this relation is of the same type as the scal- 
ing laws for the hopping regions close to the metal- 
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FIG. 1: (Color online) Two possible features of the phase 
diagram of a two-dimensional charge carrier system in the 
vicinity of the MIT: (a) At T = 0, a sample is either insulating 
(I) or metallic (M), the regions are separated by a quantum 
critical point (QCP) marked as •. However, there is a smooth 
transition between both regimes at small, finite T. (b) The 
QCP is connected with a line of continuous phase transitions 
so that the MIT should be directly detectable at finite T when 
n is varied. 



superconductor transition in ultrathin amorphous Bi 
film a 14 i 15 and close to the MIT in several disordered 
three-dimensional systems i 16 i 17 ' 18 i 19 i 20 Compare also 
Refs. Hi] and[13 for ultrathin amorphous Pd films and for 
Si-Sio.87Geo.i3-Si quantum-well structures, respectively. 

Such a scaling of the T dependences of p close to a 
phase transition, where the characteristic hopping tem- 
perature vanishes, implies the existence of a sharp bound- 
ary between apparent metallic and nonmctallic conduc- 
tion concerning p, independent of T , and hence also with 
respect to n, independent of T. Of course, this statement 
is based on some idealizations and presumes T not to ex- 
ceed a certain threshold where additional mechanisms be- 
come important. But it holds also if scaling of the T de- 
pendence is valid only in the region of activated conduc- 
tion. It will be shown below that the validity of this type 
of scaling for finite T up to the transition point would 
almost always cause a peculiarity in p(T = const., n). 

Recently, Lai and coworkers studied the transport 
properties close to the apparent MIT of an n-type Si 
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quantum well confined in a Sio.75Geo.25/Si/Sio.75Geo. 25 
heterostructure^ The authors observed the conductivity 
a to have nearly linear temperature dependences around 
the Fermi temperature Tp varying between roughly 2 and 
2.5 K, see Fig. 1 of that work. They showed that 00 (n), 
the T = conductivity obtained by linear extrapolation 
from this T region, exhibits two regimes of different slope. 
Particularly important, there is a sharp bend at the tran- 
sition between the two regimes. According to Fig. 2 of 
Ref. [H, this sharp bend could also be described as a 
knee. It coincides with the n value where da/dT changes 
its sign as T — * 0. The authors interpret this finding as 
an indication of the existence of two different phases. 

Lai et al. stress the slope of a(T) for T ~ Tf to be al- 
most constant within the n region close to the MIT. This 
statement provokes two interesting questions: (i) Pro- 
vided the extrapolated <Jo(n) has a knee indeed, and the 
slope used in extrapolation is roughly constant, should 
not a(T = const., n) exhibit a knee also for finite T? 
(ii) If yes, is the existence of this peculiarity restricted to 
the region of linear a(T, n = const.) around Tp, or is it a 
more general phenomenon? 

A first answer is already contained in Ref. I23I : In per- 
forming a percolation analysis, Lai and coworkers consid- 
ered a(T = 0.3 K, n) as best available approximation of 
<x(T = 0, n). Fig. 3 of their paper shows that deviations 
from the best percolation fitting occur just at the n value 
where oo(n) has its sharp bend. The authors considered 
this to support "that the electronic phase in the low- 
density regime is different from that in the high-density 
regime" . However, Lai et al. did not further follow up 
this line of thought but emphasized the necessity of mea- 
surements at still lower temperatures. 

The detailed consideration of the above two questions 
is the aim of the present work, being organized as follows. 
In Sec. II, reconsidering data from Ref. HH, it is shown 
that several phenomcnological aspects of a(T, n) seem to 
indicate the existence of phase transitions between ap- 
parent metallic and activated conduction at finite tem- 
peratures. A comparison with two previous experiments 
is given in Sec. III. Sec. IV is devoted to a scaling anal- 
ysis of the data from Ref. [23|. It yields an explanation 
of the shape of a(T = const., n) close to the apparent 
MIT. Finally, Sec. V summarizes and draws conclusions 
concerning the design of future experiments. 



II. PHENOMENOLOGICAL 
RECONSIDERATION OF THE EXPERIMENT 
BY LAI AND COWORKERS 

To find out to which extent and in which form the 
sharp bend of cro(n) persists when T is finite, I digitized 
Fig. 1 (b) of Ref. [23| by means of WinDig^i, considering 
the high-resolution version of this figure available from 
the cond-mat server. Values of a(T = const., n) for five 
different temperatures, 0.35, 0.64, 1.18, 2.18, and 4.0 K, 
were read out in this way. They are given in Fig. 2 
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FIG. 2: (Color online) Charge carrier concentration depen- 
dence of the conductivity for the Si quantum well structure 
studied in Ref. [23l. Data were obtained from Fig. 1(b) of that 
work. The arrow marks the critical concentration n c of the 
apparent MIT, see text. 

in a a versus n representation. First of all, this graph 
shows that all er(T = const., n) curves share a peculiar 
feature, namely an "indentation" at n « 0.32 -10 11 cm~ 2 . 
It occurs close to the concentration value where da/dT 
changes its sign at the lowest T considered in Ref. [23|, 
see Fig. 1(b) therein. This value, n c = 0.322 ■ 10 11 cm" 2 , 
is the critical concentration of the apparent MIT. It is 
marked by arrows in Figs. 2, 3, and 6 to 9 

To answer the question whether or not the sharp bend 
observed by Lai et al. at ao(n) persists at finite T, fits by 
two linear functions, concerning the n regions below and 
above 0.32 • 10 11 cm" 2 , were tried. Fig. 2 includes such 
linear fits for the case T = 4.0 K. Note, how remarkably 
well the two dashed lines describe the measured data. 

However, there is a small inconsistency in this ap- 
proximation: The two linear functions intersect at n = 
0.317 • 10 11 cm -2 , a slightly lower value than the interval 
splitting concentration. This problem is avoided when 
approximating the experimental data points by the piece- 
wise linear function f p i(n) = a + bn + c(n — n^Oin — «k) 
with the adjustable parameters a, b, c, and n^. Here, rik 
denotes the knee position, and 9 stands for the Heaviside 
step function. When included in Fig. 2, such a piece- 
wise linear fit of a(T = 4.0 K, n) would deviate from the 
dashed lines by less than the line width. So it is not 
shown for clarity reasons. 

For the piecewise linear fit of a(T — 4.0 K, n), the least 
square sum % 2 of the deviations between fit and data 
points amounts to 0.010 e 4 /h 2 when weight 1 is ascribed 
to all points. For alternative approximations by polyno- 
mials of second or third order, one obtains clearly worse 
values, x 2 ~ 0.029 e 4 /h 2 and 0.024 e 4 /h 2 , respectively. 
This is a first indication for the MIT being connected 
with a sharp continuous phase transition at finite T. 
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TABLE I: Results of phenomenological fits of a(T = const., n) 
by piecewise linear functions (plf ) and polynomials of third or- 
der (pto): Values of knee positions n k and least square sums 
X 2 are given in units of 10 11 cm -2 and e 4 /h 2 , respectively. 
The left group of columns results from consideration of the 
complete data sets, the right group is obtained from the data 
points fulfilling 0.5 a(T,n c ) < <x(T,n) < 2.0 a(T,n c ). N de- 
notes the number of data points taken into account. 
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Also the curves for the other T values can be nicely ap- 
proximated by piecewise linear functions. Tab. I presents 
the optimum knee positions ro k and the related \ 2 values. 
It is remarkable that the same n k value, 0.318 TO 11 cm' 2 , 
results for all temperatures but T = 0.35 K. This value 
is only slightly lower than the abscissa of the intersection 
for the T — > extrapolation ao(n) in Fig. 2 of Ref. l23l . 

However, when comparing with polynomial approxi- 
mations for the total n region, it is noted that the advan- 
tage of the approximation by piecewise linear functions 
over that by one polynomial of third order (same num- 
ber of adjustable parameters) diminishes with decreas- 
ing T: The approximation by piecewise linear functions 
is clearly advantageous for 4.0, 2.18, and 1.18 K, it has 
roughly the same quality as the polynomial approxima- 
tion for T = 0.64 K, and it is clearly worse than the 
polynomial approximation for T = 0.35 K, see Table I. 
This trend is not surprising because a(T = const., n) is 
basically a nonlinear function in the region of activated 
conduction, it vanishes in some exponential way with de- 
creasing n. Hence, when considering a fixed n range and 
lowering T, the a region (on a logarithmic scale) broad- 
ens, and the basic nonlincarity gets increasing weight 
compared to the "indentation" at n « 0.32 • 10 11 cm -2 . 

A more meaningful comparison of both approximations 
is obtained restricting all fits to comparable a intervals, 
namely to {0.5 cr(T, n c ), 2.0 a(T, n c )}. Under this condi- 
tion, for each of the T values considered in Fig. 2, piece- 
wise linear functions clearly approximate the data better 
than a polynomial of third order, see Table I. Moreover, 
now all fits yield the same n k value, 0.318- 10 11 cm' 2 . 

Alternative fits concerning the conductivity region 
{0.5 cr(T, n k0 ), 2.0 a(T, n k0 )} with n k0 = 0.318- 10 11 cm" 2 
yield very similar results: The advantage of the piecewise 
linear fits is confirmed for all T. The only change of the 
n k value occurs for 0.64 K where the alternative fit results 
in 0.316- 10 11 cm' 2 . 

Nevertheless, one should not overestimate the piece- 
wise linear fits. The conclusive message from their ad- 



vantage over the polynomial fits is only that, for all T 
values, some characteristic change occurs within a very 
small n region close to the knee of the piecewise linear 
fit. In this sense, the sharp bend observed by Lai et al. 
at the T — ► extrapolation <jq{ti) turns out to be not 
smoothed for finite T. Moreover, it seems to exist even 
for T < Tp, where a(T,n = const.) is nonlinear and 
considerably deviates from the estimate by means of the 
linear extrapolation from the vicinity of Tp. 

It is enlightening to compare with the behavior of 
da/dT as T -> according to Fig. 1(b) of Ref. H The 
sign change of da/dT at n c = 0.322 • 10 11 cm -2 and 
a c ~ 1.5 h/e 2 , defining the apparent MIT, roughly co- 
incides with the sharp bends in a(T = const., n). Thus 
it is likely that the bends originate from sharp continu- 
ous phase transitions separating the regions of apparent 
metallic and activated conduction at finite T. The re- 
lated critical concentration seems to be nearly indepen- 
dent of T. 

The values of n k and n c slightly differ from each other. 
The small discrepancy may arise from the approxima- 
tion of cr(T = const., n) by a piecewise linear function 
being an oversimplification. This is also suggested by in- 
specting the data points in the immediate vicinity of the 
apparent MIT in Fig. 2. 

Thus a more detailed investigation is desirable. An al- 
ternative presentation of the "indentation" with a better 
resolution is given for that purpose in Fig. 3. Here, the 
resistivity p in dependence on n is shown (so that the "in- 
dentation" is directed "upwards"). A log 10/ o scale is used 
as in many experimental papers on the apparent MIT in 
two-dimensional systems. However, this graph focuses 
on a rather small concentration region, n merely varies 
between 0.296 and 0.352 • 10 11 cm -2 . For optimum reso- 
lution, each p(T) curve is shown in a separate subgraph 
where its average slope is scaled to a common value. In 
order to make details of the curves better visible, two 
purely phenomenological fits are included: a linear re- 
gression (dashed) based on the 4 points from 0.328 to 
0.352 • 10 11 cm -2 and a regression parabola (dash-dotted) 
for the 6 points from 0.318 to 0.344 • 10 11 cm' 2 . In these 
fits, log 10 p is considered as function of n. 

Fig. 3 shows that deviations from the two fits occur 
rather suddenly, independently of T always at the same 
concentration values. Thus p(T = const., n) seems to 
exhibit three regions of differing behavior: below 0.318 ■ 
10 11 cm' 2 , between 0.318 and 0.328 • 10 11 cm' 2 , and 
above 0.328 • 10 11 cm' 2 . The medium part resembles a 
rounded step. Remarkably, its width seems to be almost 
independent of T. Its amplitude at n = 0.318-10 11 cm -2 
can be estimated by means of comparing with the linear 
extrapolation (of log w p in dependence on n) from the n 
region above 0.328 • 10 11 cm -2 . It decreases from roughly 
26 % to roughly 5 % as T increases from 0.35 to 4.0 K. 

Due to its T-independent width, the "rounded step" 
may indicate a (with respect to the experimental res- 
olution) sharp continuous phase transition occurring at 
finite T. Note that the sign change of da/dT as T — > 0, 
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FIG. 3: (Color online) Details of the dependence of the re- 
sistivity on the charge carrier concentration in the immediate 
vicinity of the MIT for Si quantum well data shown in Fig. 2. 
For the explanation of the dashed and dash-dotted lines see 
text. 

marked by the arrow, occurs roughly at the middle of 
the "rounded step" in p(T = const., n). This coincidence 
can be regarded as evidence for the hypothetical phase 
transition being the finite temperature extension of the 
zero-temperature effect metal-insulator transition. 

In order to exclude the possibility that the peculiarity 
might be pretended by the specific choice of the phe- 
nomenological regression curves included in Fig. 3, lin- 
ear regressions for the intervals from 0.282 to 0.303 and 
from 0.318 to 0.328 ■ 10 11 cm" 2 were tried as an alter- 
native. From this perspective, the a(T = const., n) data 
sets exhibit the same characteristics as in Fig. 3: There 
are three n regions, clearly differing concerning the slope 
of log 10 p(n). The apparent MIT falls in the middle re- 
gion, the width of which is roughly constant again. This 
supports the above proposed interpretation in terms of a 
phase transition between activated and apparent metallic 
conduction at finite T. Because of its similarity to Fig. 
3, the corresponding graph is not shown here. 



III. COMPARISON WITH OTHER 
EXPERIMENTS 

One could object the above analysis furthermore with 
the following argument: Even if the peculiarity under 
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FIG. 4: (Color online) log-log diagram of p(T = const., n) 
for the intermediate vicinity of the apparent MIT in a high 
mobility MOSFET, reproduced from Fig. 2 of Ref. [H Data 
were obtained from Fig. 1 of Ref. HH , The dashed lines serve 
only as guide to the eye. 
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FIG. 5: (Color online) log-log diagram of p(T = const., n) for 
the vicinity of the apparent MIT in an AlAs quantum well, 
reproduced from Fig. 3 of Ref. [25j. Data were obtained from 
Fig. 2 of Ref. [2^, considering the high-mobility direction. The 
dashed lines are included as guide to the eye. 

discussion is clearly visible at five different temperatures, 
it might be an artifact of some systematic experimental 
imperfections specific to Ref. [23[ This objection could 
also be formulated as the question whether or not similar 
peculiarities have been reported in previous experimental 
studies. In fact, to the best of my knowledge, this is not 
the case. 

However, a reconsideration comparing experimental 
p(T) data for the apparent MIT in two-dimensional sys- 
tems with experience on three-dimensional disordered 
systems was performed five years ago.— It points to possi- 
ble peculiarities in p(T = const., n) close to the apparent 
MIT in measurements at an high mobility Si MOSFET^ 
by Kravchcnko et al. and at an AlAs quantum well 2 ^ by 
Papadakis and Shayegan. The respective graphs from 
Ref. are reproduced here as Figs. 4 and 5. 
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Fig. 4 shows that already in the results of Ref. [HI, 
one of the first studies of the apparent MIT, a behav- 
ior may be hidden which is similar to that in Fig. 3: 
The p(T = const., n) curves may exhibit a shoulder-like 
structure close to p = h/e 2 being related to an offset 
of the curve parts below and above h/e 2 . This struc- 
ture resembles the "rounded step" in Fig. 3 discussed 
above. It is located close to the intersection region of the 
p(T = const., n) curves, for a more detailed discussion 
see Ref. H^. Also the AlAs quantum well studied by Pa- 
padakis and Shayegan^ seems to exhibit such an offset 
in p(T = const., n) close to the apparent MIT as Fig. 5 
displays. 

One point is particularly remarkable: For the lowest 
considered T, the highest resistivity values of the high n 
branch of p(T = const., n) are roughly the same in the 
three experiments considered here. They amount to 0.78, 
0.80, and 0.63 h/e 2 for Refs-HQJ, andU, respectively. 

Concluding the comparison, Figs. 4 and 5 resemble Fig. 
3 substantially, but they are certainly not as conclusive 
as the reconsideration of data from Ref. above. Nev- 
ertheless, in my opinion, they give additional support to 
the hypothesis of sharp continuous phase transitions at 
finite temperatures. 

For all the three experiments reconsidered here, one 
could of course object that the emphasized features might 
arise from random deviations in the measured values of 
n and p. However, it seems to be very unlikely that all 
the following coincidences occur only by chance: 

(i) These features seem to be present for all tempera- 
tures considered. 

(ii) They are observed at data from three independent 
experiments with different kinds of samples, made up of 
different materials. 

(iii) The peculiarities have qualitatively the same form 
in all three data sets. 

(iv) They occur in all cases in the same resistivity re- 
gion, slightly below h/e 2 . 

(v) They are observed in all cases close to the common 
intersection point of the p(T = const., n) curves. 

It has to be mentioned that similar hints to a possi- 
ble peculiarity could not be found in data from other 
publications, for example Ref. However, this is not 
a strong objection to the reconsideration presented here 
for the following reasons. 

Possibly, on the one hand, very low charge carrier den- 
sities together with extremely high mobilities as realized 
in Ref. [23| might be needed. On the other hand, the pe- 
culiarities may be easily overlooked: Many of the exper- 
imental studies in the literature concern broad n ranges, 
and do not obtain so many data points in the immediate 
vicinity of the apparent MIT, at both sides, as Ref. |23[ 
Thus the peculiarity is small on the logarithmic p scale 
usually used and may easily be covered by the nonlinear 
n dependence of a, compare the effect of restricting the a 
range discussed in Sec. II. Moreover, when considering p 
on a logarithmic scale, there is no need to achieve a very 
high precision of the individual data points, so that the 



small peculiarity may remain hidden behind random de- 
viations. Additionally, arbitrary sample inhomogeneities 
including edge effects tend to smear it, see the careful 
study in Ref. [28|, in particular Fig. 4. The last three 
reasons seem to me more likely explanations than the 
possible necessity of an extremely high mobility. 

IV. SCALING ANALYSIS 

Provided the "rounded step" is indeed related to a 
phase transition, how may it arise? On a macroscopic 
level, it can be explained in the following way: Suppose, 
on the insulating side of the apparent MIT, the T depen- 
dence of a scales as observed by Kravchenko et al. at the 
resistivity of MOSFET's^ 

a(T,n) = a(t) with t = T/T (n) , (1) 

for further references concerning this matter see the in- 
troduction. To the best of my knowledge, there is no sys- 
tematic study how this scaling is softened with increasing 
T close to the MIT. So, according to the above hypothe- 
sis, assume that it holds up to the critical concentration 
n c . Since da/dT — > as n — > n c , 2o(n) vanishes and t 
diverges as n — > n c , independent of T. This is connected 
with cr(T, n c ) = a(t = oo) = er c (provided T does not ex- 
ceed a certain threshold where other mechanisms become 
relevant). 

Suppose furthermore 

T (n) = A\6nf with 5n = n-n c , (2) 

where A is a constant which might, however, be sam- 
ple dependent. The value of the critical exponent was 
determined by Kravchenko et al. by means of a scal- 
ing analysis of their MOSFET experiment yielding (3 = 
1.6 ± O.lM- This result was confirmed by Lam et al. for 
a Si-Sio.87Ge .i3-Si quantum well, j3 = 1.6 ± 0.2^ (The 
comparison with the behavior of ultrathin metal films an- 
alyzed in Ref. [TH is hindered by the use of the resistance 
as control parameter in that work.) 

According to Eqs. (1) and (2), a should depend only 
on the parameter T/\5n\^ or rather on dn/T 1 ^. This 
can be checked analogously to e.g. Fig. 10 of Ref. [H or 
Fig. 3 of Ref. [H. Since the validity of Eq. (2) is needed, 
such a scaling check is less general than the construction 
of a mastercurve by rescaling of T. However, it has a big 
advantage: Natural small sample inhomogeneities An do 
not destroy a possible generic coincidence of the curves 
for \Sn\ 3> An, although such inhomogeneities might 
even cause nonmonotonic behavior of a(T,n = const.) 
in the immediate vicinity of n c . 

A check of the hypothesis that a is a function of 
on /T 1 /' 9 alone is given in Fig. 6 considering the a curves 
for T = 0.35 and 0.64 K from Fig. 2. For that plot, 
n c = 0.322 • 10 11 cm~ 2 corresponding to the sign change 
of da/dT as T — > 0, and /3 = 1.6 as obtained in Ref. [l3| 
arc utilized. Curves for higher T are not taken into ac- 
count since cr(T, n c ) becomes T dependent above roughly 
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FIG. 6: (Color online) Check regarding scaling behavior ac- 
cording to Eqs. (1) and (2) for the two a curves from Fig. 
2, related to the lowest T values, 0.35 and 0.64 K. Note: 
The agreement of both curves does not result from adjusting 
parameters. 
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FIG. 8: (Color online) Scaling check with respect to Eqs. (2) 
and (4) for all a curves presented in Fig. 2, presuming j3 = 2.0. 

phous Sii-^Cr^ filmsi^ that, close to n c , the multiplica- 
tive decomposition 
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FIG. 7: (Color online) Scaling check with respect to Eq. (2) 
and (4) for all a curves presented in Fig. 2, presuming p — 1.6. 



1 K, presumably because some additional mechanism be- 
comes relevant there, see Fig. 1 of Ref. [23|. 

Fig. 6 shows that both a curves nicely fall together. 
Note also the conformity concerning the "indentation" 
around Sn = 0. Since the transformation used here does 
not include any adjustable parameter, it is unlikely that 
the agreement of both curves arises only by chance. Thus 
Fig. 6 supports the scaling equations (1) and (2). 

Now the question is suggested whether signs of the 
scaling demonstrated in Fig. 6 can be detected also 
for T > 1 K where a(T,n c ) has a significant T de- 
pendence. Here, a comparison with the MIT in three- 
dimensional systems is helpful. It was observed at amor- 



a{T, n) = a sca i(T/T (n)) ■ £(T, n) (3) 

is useful, where <7 sca i denotes the scaling contribution 
valid in the low-temperature region, and £, with £ — > 1 as 
T — > 0, describes the hi gh-T deviations from the scaling 
behavior. Fig. 5 of Ref. |l8| shows that £(T, n) proved to 
be almost independent of n in the vicinity of the MIT. 
Presuming for the moment that £ is indeed independent 
of ri, Eq. (3) implies 

a{T,n)/a(T,n c ) = cr sca ,i(T/T (n))/cr c . (4) 

In this sense, for the data considered here, the ratio 
a(T,n)/a(T,n c ) might be a function of dn/T 1 ^ alone 
within the whole T region from 0.35 to 4.0 K. This 
hypothesis is checked in Figs. 7 and 8 assuming n c = 
0.322 • 10 11 cm" 2 as in Fig. 6. 

Fig. 7 is based on j3 = 1.6 as Fig. 6. It shows in- 
deed an approximate collapse of the data for various T 
on one curve: Except for T = 4.0 K, all points with 
Sn/T°- e25 < 0.03 • 10 11 cm- 2 K-°' 625 can be considered 
to fall onto one curve. However, for larger values of 
Sn/T 625 , there are clear deviations. They can have sev- 
eral causes: (a) The scaling equation (1) is expected to be 
violated somewhere in the region of metallic conduction. 
Otherwise, a(T = 0, n) would have to be independent 
of n. (b) The deviations may originate from a weak n 
dependence of £. (c) Inaccuracies of the values of n c and 
(3 may be the reason. 

Fig. 8 demonstrates that a shift of (3 from 1.6 to 2.0 
leads to a far better data collapse. (Its quality is once 
more considerably improved if the points for the lowest 
temperature, T = 0.35 K, are omitted.) However, it re- 
mains an open question at present whether the improved 
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agreement of the curves has a physical meaning or arises 
from the compensation of several deviations from scaling. 

Nevertheless, in collapsing points for various T values 
onto one curve, Figs. 6 to 8 clearly testify that, inde- 
pendently of T, all a(T = const., n) curves exhibit a 
sharp bend very close to n c . Thus they support once 
more the hypothesis of a line of phase transitions at fi- 
nite T. Simultaneously, the above consideration has the 
consequence that cr(T, n c ) has the meaning of a minimum 
metallic conductivity also for T > 1 K even though it is 
T dependent in this region. 

Although the argumentation in this chapter is solely 
based on scaling according to Eq. (1) within the region 
of activated conduction, the behavior on the apparent 
metallic side as shown in Figs. 6 to 8 likewise deserves 
attention. These graphs, in particular Fig. 6, may indi- 
cate analogous scaling behavior also for a certain n re- 
gion above n c , in agreement with the original findings in 
Ref. Il~3l . but in contradiction to later publications. Since 
such a scaling would not be consistent with usual metal- 
lic transport, systematic experimental studies of how this 
scaling breaks down should be very interesting. 

The quantitative determination of the critical expo- 
nent (3 is hampered by the lack of a theoretically based 
analytic expression for cr sca x(t) being applicable to a 
broad a region. The only resort is to use a sufficiently 
precise empirical description of the hypothetical master 
curve fiSn/T 1 ^) = r(T,Sn) = a(T,n)/a(T,n c ) with an 
as small as possible number of adjustable parameters. 
In practice, for several temperature pairs (Ti,T2), the 
value of the scaling factor X = (Ti/T^) 1 /' 3 is obtained 
by "moving" the experimental curve r(T2, X5n) on top 
of the curve r(Ti,Sn). This is realized approximating 
r(Ti,5n) and r(T2, XSn) by a common empirical g(Sn) 
with several adjustable parameters. In minimizing the 
total mean square deviation of both measured data sets 
from g(6n), the parameters of g and the scaling factor X 
are adjusted. 

This procedure is complicated by the sharp bend in 
er(T = const., Sn): For example, a polynomial would have 
to be of a very high degree to emulate this feature. A 
way out is to "uncompress" the immediate vicinity of n c 
by considering z = sign(<5n,)|i5n| 1 / 2 as independent vari- 
able instead of Sn. Thus, as various trials corroborated, 
r(T, z) = r(T, Sn) can be well approximated by a g(z) be- 
ing a polynomial of rather low degree or a natural cubic 
spline with a rather small number of equidistant nodes. 

Table II presents values of the scaling factor X ob- 
tained in a series of such polynomial regression analy- 
ses. The errors of X given are rough estimates of the 
lcr bound presuming that the deviations from g(z) are 
only of random nature, and that the relative errors have 
Gaussian distributions of same width for all data points 
taken into account. For four pairs (Ti,T2), the table 
compares the result of investigating the total r(T, z) data 
sets with the two values obtained when omitting the four 
outer points with |<5n| > 0.04 ■ 10 11 cm" 2 from each ex- 
perimental curve, and when considering only the region 



TABLE II: Results of scaling regression analyses determining 
the best overlap of the curves a(Ti,n c + Sn)/a(Ti,n c ) and 
a(T2,n c + X 5n)/a(T2, n c ). The three columns of X values 
refer to the study of complete a(T, n) data sets, to the analysis 
of data sets restricted by the demand \Sn\ < 0.04 ■ 10 11 cm -2 , 
and to the consideration of only the points within the region 
of activated conduction, \5n\ < 0. 



Ti(K) 


T 2 (K) 


X(compl. set) 


X(restr. set) 


^(activ. cond.) 


0.35 


0.64 


0.684 ± 0.017 


0.693 ± 0.027 


0.646 ± 0.037 


0.64 


1.18 


0.719 ± 0.015 


0.729 ± 0.021 


0.686 ± 0.045 


1.18 


2.18 


0.744 ± 0.011 


0.759 ± 0.016 


0.770 ± 0.037 


2.18 


4.0 


0.788 ± 0.012 


0.811 ± 0.016 


0.781 ± 0.038 



of activated conduction, \Sn\ < 0. In these three cases, 
polynomials of degree 5, 4, and 3, respectively, were used 
as g(z). Additional analyses based on polynomials of by 
1 to 2 enlarged degree yielded very similar results, as did 
studies based on natural cubic splines. 

Table II shows that, for all (Ti, T2) pairs, the three X 
values agree nicely with each other. Since the r(T,5n) 
curves have a considerable curvature and, moreover, ex- 
hibit the sharp bend fine structure, this agreement sup- 
ports the phcnomcnological Eqs. (1) to (4) within exper- 
imental accuracy in the following sense: The n depen- 
dence of £ seems to be very weak indeed, otherwise the 
values in the first two X columns would have to differ sys- 
tematically. Moreover, the agreement between first and 
third X columns supports the hypothesis obtained above 
from Figs. 6 to 8 that scaling behavior may be valid also 
in a certain n region of apparent metallic conduction. 

The question remains whether the deviations from 
scaling in Fig. 7 may arise mainly from using an im- 
precise value of /3, or whether Figs. 7 and 8 point to 
systematic deviations from scaling indicated by a feigned 
drift of (3. This exponent can be directly obtained from 
the X values in Table II. However, in estimating its error 
bound, one has to have in mind that the measurements in 
Ref. l23j were performed in a dynamic modus with slowly 
varying temperature, sec Fig. 1 therein. Additionally to 
calibration uncertainties, such a procedure exhibits some 
small systematic error of the T values arising from im- 
perfect equilibration. Even though these inaccuracies are 
so small, that they are extraneous to the data analysis 
in Ref. |23| . they may have a significant influence on the 
(3 values obtained here. Thus, to judge the influence of 
systematic T errors, the corresponding uncertainties of j3 
were estimated, cautiously guessing AT/T ~ 0.03. 

In this way, the /3 values 1.6 ± 0.5, 1.9 ± 0.5, 2.1 ± 0.5, 
and 2.5 ± 0.7 are obtained considering the intervals 
0.35-0.64 K, 0.64-1.18 K, 1.18-2.18 K, and 2.18-4.0 K, 
respectively. Here, the error originates from the 3cr un- 
certainty of X and from the systematic T deviations. Ac- 
cording to these results, the effective (3 value might slowly 
vary with T due to some deviations from the very simple 
Eqs. (1-4). However, for the size of the error bars, this 
variation of (3 can still not be taken for certain. Further 
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FIG. 9: (Color online) d log 10 a /An versus n for all curves 
displayed in Fig. 2. 

experiments are needed to clarify this question. 

We focus now again on the region of activated conduc- 
tion. For the immediate vicinity of the MIT, that means 
for large t and \a c — a\/a c -C 1 , the following ansatz is 
suggested: 

o{t) = a c -(l-Br v ) (5) 

Here B is a dimensionless constant and v a positive ex- 
ponent. Unfortunately, only a very rough guess of the 
value of v can be given at present. If the value would 
correspond to hopping in the Coulomb gap for small t, 
v would be 1/2. However, according to Ref. [13, v may 
be expected to be smaller for large t. In both cases, for 
/3 = 1.6 according to Ref. [TH, the product j3v would be 
clearly smaller than 1. 

Equations (1), (2), and (5) determine the type of the 
n dependence of a close to a c : 

a(T,n)=a c -{l-C\6nf u ) (6) 

where 

C = B{A/T) V . (7) 

Thus, only if (3v = 1, a{T = const., n) does not exhibit 
any specific feature when n approaches n c from the side of 
activated conduction; nevertheless, a(T = const., n) may 
have a knee at n c in this special situation. In case /3v < 1 
as guessed above, a(T = const., n) has a root-like pecu- 
liarity at n c , and p(T = const., n) as well. Taking into 
account small sample inhomogeneities, this peculiarity is 
softened to some "rounded step" . This is in qualitative 
agreement with Fig. 3. 

A root-like peculiarity in <j{T = const., n) would im- 
ply a divergence of dcr/dn, and thus also of dlog 10 <r/dn. 
(Considering log 10 a rather than a has the advantage that 
the contribution from £, see Eq. (3), is suppressed.) The 



divergence would be smoothed to a sharp maximum in 
reality, caused by sample inhomogeneities and by trunca- 
tion errors in the numerical evaluation of o(T = const., n) 
utilizing a difference quotient. 

The hypothesis of the sharp maximum in dlog 10 <7/dn 
is checked here using the same approach as in Eqs. (3), 
(8), and (9) of Ref. [30: A window including 4 neighboring 
points is shifted along the log 10 a curve. The slope values 
from the corresponding linear fits of log 10 a as function 
of n are related to the n values, for which polynomials of 
second order would be exactly differentiated. To consider 
4 points in the fits turned out to be a good compromise 
between the demands for high resolution and for low ran- 
dom errors. 

Fig. 9 shows that, independent of T, dlog 10 (r/dn ex- 
hibits the expected sharp maximum - note that a loga- 
rithmic scale is used to display the derivative values -. 
It is striking that the peak is always located just at n c . 
This feature is another indication of phase transitions 
occurring at finite T when n crosses n c . 

From the phenomenological point of view, the peak of 
the derivative has two effects: On the one hand, it is the 
reason of the small differences between n c and the rik 
values from the piecewise linear fits, see Sec. II. On the 
other hand, it makes the sharp bend in er(T = const., n) 
more easily visible. 

The question remains, whether the T dependence of 
the amplitude of the rounded step can also be un- 
derstood. As discussed above, the amplitude of the 
"rounded step" decreases with increasing T as expected 
from Eqs. (6) and (7). However, for a quantitative anal- 
ysis, more detailed experimental data are needed. 

Finally, a counterintuitive consequence of the scaling 
is worth pointing to: In measuring <r(t), the influence of 
sample inhomogeneities decreases with increasing T since 
the n range needed to explore a certain fixed a interval 
broadens. Hence measurements at intermediate T should 
be more promising than studies at extremely low T . 



V. CONCLUSIONS 

In summary, reconsidering the experiment of Lai and 
coworkers,— a series of features has been described which 
indicate that the apparent MIT at T = may be con- 
nected with a line of sharp continuous phase transitions 
at finite T: (a) Close to the critical charge carrier concen- 
tration n c , defined by the sign change of da/dT as T — > 0, 
piecewise linear functions approximate <j(T = const., n) 
clearly better than polynomials of third order. This holds 
for the entire T range from 0.35 to 4.0 K. (b) The knee 
of the piecewise linear function is always located close to 
n c . (c) In the immediate vicinity of n c , p(T = const., n) 
has a "rounded step" structure of T-independent width 
what is connected with an offset in the curves, (d) The 
middle of the "rounded step" always coincides with n c . 
(e) Two previous experiment s 13 ! 26 exhibit similar offsets 
in p(T = const., n) close to n c . (f) The offsets in the 



9 



data from Rcfs. [l3ll23l . andl26l occur at roughly the same 
p value, a little below h/e 2 . (g) The data from Ref. I23I 
exhibit sharp peaks in the derivative d log 10 cr/drt at n c 
for all considered T values. 

Moreover, it has been shown that the "rounded step" 
in the p(T = const., n) curves can be explained in terms 
of scaling of the T dependences of a for various n within 
the region of activated conduction. The applicability of 
this scaling is supported by the collapse of the curves 
for T = 0.35 and 0.64 K when plotting a as function of 
(n — n c )/T 625 , utilizing the value of the critical expo- 
nent of the characteristic temperature obtained in Ref. 
[iH Additionally, an approximate collapse is observed for 
the whole T region from 0.35 to 4.0 K when considering 
the ratio a(T,n)/a(T,n c ) rather than a(T,n). 

The scaling of cr(T, n) for the region of activated con- 
duction should cease to be valid when the characteristic 
hopping temperature reaches zero indicating a qualita- 
tive change in the transport mechanism. This happens 
at the T-independent resistivity value, a c = a(T,n c ). In 
this sense, scaling gives additional, indirect support to 
the hypothesis of a line of finite temperature phase tran- 
sitions from activated to apparent metallic conduction. 

As a whole, the above indications seem to provide 
rather convincing evidence for a line of sharp phase tran- 
sitions at finite T. Nevertheless, to the best of my knowl- 
edge, there is no microscopic theory available explaining 



this feature. So further theoretical work as well as addi- 
tional experiments in the immediate vicinity of the ap- 
parent MIT are called for. These experiments should pri- 
marily focus on enhancing precision, in particular with 
respect to homogeneity of the samples, rather than on 
reducing the lowest accessible temperature. If they turn 
out to support the above conclusion, they will shed new 
light on the problem whether or not the apparent MIT is 
a real phase transition, and set qualitative constraints on 
theoretical models. In such investigations, scaling anal- 
yses as performed above may be helpful to clarify the 
character of the apparent metallic state. 
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